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Abstract 

We demonstrate that for a migration of a liquid layer between the melting and 
the solidification front an exact steady-state solution with two parabolic fronts can 
be found. A necessary condition hereby is that the temperature of the solidification 
front exceeds the temperature of the melting front (both temperatures are supposed 
to be constant). It is shown that in pure materials and alloys there exist two types 
of solutions with two convex and with two concave parabolas respectively. While a 
steady-state process with two planar interfaces is only possible for a single point, the 
processes with two parabolas are possible inside a region of control parameters. The 
relations between the Peclet numbers and the control parameters are obtained. 
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1 Introduction 



Liquid film migration (LFM) is well-known phenomenon which has been observed in many 
alloy systems during sintering in the presence of liquid phase an< ^ m Cu-In solid solutions 
during melting started at grain boundaries jS]. In LFM, one of crystals is melted and the 
other one is solidified. The both solid-liquid interfaces are moving together with the same 
velocity. In the investigated alloys systems the migration velocity is of the order of 10~ 8 — 10~ 7 
m ■ s _1 and it is controlled by the solute diffusion through a thin liquid layer between the 
two interfaces j3] . The migration velocity is much smaller than the characteristic velocity of 
atomic kinetics at the interfaces. Therefore the both solids at the interfaces must be at the 
local thermodynamic equilibrium with liquid. The theory [3| answers the question about the 
different equilibrium states at the melting and solidification fronts. In a steady-state regime, 
the difference in equilibrium relates to coherency stresses appearing only at the melting front 
due to the sharp profile of composition ahead the moving melting front. Thus, the liquid 
composition at the melting front which depends on the coherency strains and a curvature 
of the front differs from the liquid composition at unstressed and curved solidification front. 
The migration velocity is proportional to the difference of these compositions divided by the 
film thickness jlj. But what controls the thickness? 

Consequently, a new problem of two combined moving solid-liquid interfaces with a liquid 
between them appears. In the present article, this problem is considered under simplified 
boundary conditions: the temperature and chemical composition along each interface are 
constant. These constants are different for the melting and solidification fronts and differ 
from those far from the migrating liquid film. It means that any capillary, kinetic and 
crystallographic effects at the interfaces are neglected. It is found that under these simplified 
boundary conditions two co-focal parabolic fronts can move together with the same velocity. 
The situation is rather similar to a steady-state motion of one parabolic solidification front 
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into a supercooled melt |H| or one parabolic melting front into a superheated solid. 



2 Solutions for one parabolic front 

Needle-like stationary solutions were first obtained by Ivantsov for the crystallization of a 
pure material from a supercooled melt j3] and were extended to binary alloys [Hj- This 
solutions describes a parabolic interface of a solid phase at constant temperature T m which 
extends into a supercooled melt. Inside the melt the temperature drops and reaches far from 
the interface its asymptotic value T < T m . 

For the two-dimensional case the Ivantsov relation is 

A = F(P) = v / ^Pe p erfc(v / P) , < A < 1. (1) 

It connects the Peclet number P = VR/2D and the supercooling A = (T m — T )c/q. Here 
V and R are the front-velocity and the tip-radius of the parabola respectively; D and c are 
the thermal diffusivity and the specific heat of the melt, respectively, q is the latent heat. 
Eq. (^) can easily be obtained in 2D-parabolic coordinates (77, £) (see, for example, Ref. 

0) 

77 = Vx 2 + z' 2 + z , £ = \/x 2 + z' 2 - z , 0<(?7,£)<oo, 

x = ±\fr^ , z'^-{r]-^) , -00 < (x, z) < 00, (2) 

where (x, z' = z—Vt, t) are the Cartesian coordinates and the t is the time. In this description 
the parabolic solidification front [z' — ^ — |^), moving with velocity V in the positive z- 
direction, has the coordinate t] = R. The regions < 77 < R and R < 77 < 00 correspond to 
the solid and the liquid phase respectively. 

In this set of coordinates, the temperature field T in both phases depends only on 77. The 
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thermal diffusion equation 



is easily solved by the ansatz 

T(7?) = A + B J r]- 1/2 e- V7l/2D d7] (4) 

with different constants A and 5 for the solid and the liquid phase. These constants and Eq. 
(0) can be obtained by using the appropriate boundary conditions, namely the asymptotic 
boundary condition T(rj — > oo) = T and the interfacial conditions of temperature continuity 
T(R — 0) = T(R + 0) = T m , and heat balance at the interface 

Vq/c= -2DT'(R + 0) + 2DT'(R-0), (5) 

where the prime denotes the derivative with respect to 77. For reasons of simplicity we 
assumed the same values for D and c in both phases. 

Apart from the solution with convex parabolic front there exists as well a solution with 
a concave front (z' = —^ + |^). In the parabolic coordinates, the interface is defined by the 
relation £ = R while solid and liquid phase are at 00 < £ < R and R > £ > 0, respectively. 
In contrast to the former case the temperature field T depends only on £. Instead of Eq. (J3J) 
the diffusion equation therefore reduces to 

°(j- ( %)v(j- ( %) =0 , (6) 



and is solved by the ansatz 

T(0 = A + B J C 1/2 e n/2D d£ . (7) 

While the conditions of continuity and heat balance at the interface stay the same as for the 
convex parabola, the asymptotic condition changes to T(£ —►())= Tq. Thus for the concave 



front instead of Eq. (JI} we obtain the relation 

A = $(P) = 2VP e~ p e x dx < A < 1.284... , (8) 

Jo 

between the supercooling A and the Peclet number P. The function $(P), in contrast to 
the function F(P) in Eq. (JTJ, is not monotonous: it has a maximum 1.284... at P = 2.25... 
and approaches 1 from above for P — >• oo. 

Naturally, Eq. (0) and Eq. (JHJ) describe as well a parabolic melting front: the superheated 
solid phase with T > T m is outside the convex parabola (rj = R) at R < rj < oo or inside 
the concave one (£ = i?) at i? > £ > 0. The melting heat is — g. 

3 Two parabolic fronts in a pure material 

It turns out, that an exact steady-state solution with a combined motion of the two co- focal 
parabolic fronts (melting and solidification front with a liquid layer in between) can be found. 
For this the solidification front temperature T m has to exceed the temperature T m on the 
melting front. This may be the case for pure materials where the melted solid phase (SI) has 
with respect to the solidified solid (S2) an additional contribution e to the energy- density 
(for example due to preliminary mechanical treatment). In this case the melting heat equal 
to — (q — e) and the melting temperature T m = T m (l — e/q) differ from the values q and T m 
for the solidification front. The additional energy e is the driving force for a recrystallization 
process, which takes place in the solid either due to a migration of a boundary between a 
SI- and a S2-grain or due to a possible recrystallization through a liquid layer between SI 
and S2. 

It should be pointed out, that a combined stationary motion of two planar interfaces (with 
equilibrium temperatures T m and T m ) is only possible at the temperature T 0j o at which the 
temperature increase (T m — T 0j o) exactly compensates the additional energy, c(T m — T 0j o) = e. 
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All temperatures deviating from TJ^o lead to an increase (at To > To o) or to a decrease and 
disappearance (at T < T 0j o) of the liquid layer between the planar interfaces. Here only two 
co-focal parabolic interfaces can lead to a steady-state motion. 

Similar to the considerations described for one interface an analysis of the temperature 
field in parabolic coordinates leads to the following relations: 
For two convex parabolas with P\ > P 2 (Fig. 1): 



A, 



A 



2JP 2 e P2 / e~ x 



dx 



IP, 



1 - a A m - \l — e 
1 1 



2 P2~Pl 



For two concave parabolas with Pi < P2 (Fig. 1] 



(9) 
(10) 



A r 



A 



1 — p rv** 2 

2JP 2 e- P2 / e x dx, 

1 - aA m - W — e - 
V "1 



(11) 
(12) 



with Pi = VRi/2D, A m = (T m - f m )c/g, A = (T - f m )c/q and F(P) and $(P) as in 
Eq. (HJ) and Eq.(jHJ). a = q/cT m is a material parameter (e.g., for Ni one has a = 0,25). 
The term (1 — aA m ) arises due to the fact that the melting heat of a solid phase SI with 
an additional energy density e differs from the equilibrium melting heat, — q (note that 
-(<? - e ) = -Q (1 - «A m )). 

In the corresponding limiting cases Eqs. (j9j)- (fT2~j) reduce to Eqs. (JTJ) and (JHJ): Eq. Q 
for P\ — > 00 and Eq. (fTTj) for Pi — ► lead to the solidification relations Eq. (0) and Eq. 
(JHJ), respectively, with supercooling A m instead of A. Eq. (fTUj) for P 2 — > and Eq. (fT2"j) for 
P2 — ^ 00 lead to the melting relations Eq. (JTJ and Eq. (jHJ) with the normalized superheating 
A/ (1 — aA m ) instead of A. 
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4 Analysis of convex and concave solutions 



Even if Eqs. (JHJ) - (jEj) are valid for arbitrary values of the normalized driving force A m > 0, 
we consider in the following analysis the case of small driving forces A m <C 1. With this 
assumption Eqs. and (fTTj) . respectively, simplify to 



1 

2 
1 

2 



l P 1 ± y/Px - 2A r 

f P 1 + Jp 1 + 2A r 



(Pi > P 2 ) 
(Pi < P 2 ) 



(13) 
(14) 



Eq. ()13j) together with Eq. (fTU|) and Eq. ()14j) together with Eq. f!12|) define for the case of 
small driving forces A m <C 1 the dependency of A (Pi) for both types of solutions (concave 
and convex parabolas) (see Fig. 2). 

For Pi 3> A m the upper branch of Eq. (fT3|) gives 



P 2 Pi - A m (P x > P 2 ) 



(15) 



and from Eq. (|14jl 

P 2 ^Pi + A m (Pi<P 2 ). (16) 

Using this simplification and taking the limit A m — > we obtain from Eqs. (|10|) and (|12|) 
the simplified expressions 



A, 
A 



1 - a + 
1 — a — 



1 



2P 
1 



F(Pi) , 



(Pi > p 2 ) 
(Pi < p 2 ) 



;i7) 



A m V 2 Pi . 

A stationary solution with two convex parabolas (Pi > P 2 ) exists for 



with 



An < A < A* 



An = f 1 - a) A m and A* = 1 - aA r 



(19) 



(20) 



The point A = A corresponds to the steady-state solution of two moving planar interfaces 
mentioned above. 

If A < A < Ai, the dependence of A on P 1 is defined by Eq. (HTJ) and P 2 = Pi - A m . 



At the point A = A x = yVA m /2 one has Pi = 2A m and P 2 = A m /2. In the vicinity of the 
Appoint one can obtain 

s T J^~l. (21) 
From Eqs. and f|l Ojl . one can obtain that close to the limiting point A* one gets 

Pl - 2(A^A) >y 1 aml fi " ^ <<L (22) 

These relations describe the independent motion of both parabolic fronts propagating with 
the same velocity. 

Solutions of Eqs. (fT2*j) - (fUjl for a pair of concave parabolas exist in the region 

- I A 4 | < A < A 2 (23) 
between two extreme points on the A(Pi)-dependence (Fig. 2): minimum point | A4 | = A m 



at Pi = v/A m /4(l — a) and maximum point A2 which is close to Ao- In the limit Pi — >• 00 
the value of A tends to the same point Aq as for the convex parabolas, A = Aq — a/2P\. 



5 Parabolic solutions for binary alloys 

Extending our analysis to a two-component alloy, we have to take into account that in 
contrast to a pure material the temperatures T\ and T 2 of the melting and solidification 
fronts are unknown (in a pure material T\ = T m , T 2 = T m ). In order to define T\ and T 2 
(and consequently A m = (T 2 — Ti)c/q, and A = (T — T^c/q in Eqs. (pi-([12j)) we have to 
obtain two additional equations. 
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In our description, we denote the molar fraction of the second component with C and the 
diffusion coefficients in the solid and the liquid phase with D$ and Dl respectively. While 
for convex parabolas the concentration fields, C(r),£), only depend on 77, they depend only 
on £ in the concave case. In the first case, the field C(rj) satisfies the equilibrium boundary 
conditions at both interfaces, namely the continuity condition 

C(R 2 - 0) = C S (T 2 ) , C(R 2 + 0) = C L (T 2 ) , C(R 1 - 0) = C L {T X ) , C{R X + 0) = C s (Ti) , 

(24) 

the far field condition C(rj— >oo) = Co and the conservation conditions at the interfaces 



V[C L {T 2 )-C S {T 2 )} 



V 



-2D L C'(R 2 + 0) + 2D S C'(R 2 - 0) 
-2D L C\R 1 - 0) + 2D s C\R l + 0) 



(25) 
(26) 



Here Cs(T 2 ) and Cl(T 2 ) are the solidus and liquidus compositions which are defined 
by the equilibrium phase diagram of the alloy, while Cs(Ti), Cx(Ti) are the corresponding 
values defined by a disturbed phase diagram. The additional energy density e, changing 
the melting point of the melted solid phase SI, will change the equilibrium compositions 
at the melting interface as well. In this case, the composition differences Cs(T\) — Cs{Ti), 
Cl{T{) — Ci(Ti), are proportional to e/q. Another reason for the distortion of the phase 
diagram are coherency stresses due to compositional inhomogeneities in front of the melting 
interface jl]. 

By solving the diffusion equation for C(rj) and applying the appropriate boundary con- 
ditions ()24 )) -([26 [) . we obtain for the case of convex parabolas the solution: 



C L (T 2 ) - CLjTC 
C L (T 2 ) - C S (T 2 ) 

C L (T 2 ) - C S (T 2 ) 



2\jP 2L e p ^ 

ClPV - C 8 (Ti) 
C L (T 2 ) - C S (T 2 ) 
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e x dx 



P 2L p 2L ~p 1L 

PlL 



PlL > P2L 

F(Pis) 



(27) 
(28) 



(P iL = DPi/D L , Pis = D L P 1L /D S ; F(P) is the same function as in Eq. (JTJ)). A similar 
consideration for the case of concave parabolas with concentration field C(£) leads to: 

C L (T 2 ) - C L {T X ) , - rVPE 



C L (T 2 ) - C S (T 2 ) 
C L (T 2 ) - C S (T 2 ] 



2JP 2L e 



-P2L 



e x dx 



2L p PlL-P2L 
\l 



PlL < PlL 

®(Pis), 



(29) 
(30) 



C L (T 2 )-C S (T 2 ) 

where 3>(P) is defined by the same manner as in Eq. (JSJ). 

For the convex configuration, Eqs. ©-((TDJ) together with Eqs. (J2*7|) -(|2*5 |) describe the 
dependency of the quantities T\, T 2 , Pi, P 2 on the initial control parameters T , Co, and 
e. The corresponding relations for the concave configuration are given by Eqs. (|11 |) -([12 )) 
and Eqs. (|29) ) -(|30J ) . In systems with vanishing e the coherency strain effect can support the 
combined motion of two interfaces as, for example, in liquid film migration 



6 Discussion 

The considered two fronts process (TFP) is an alternative to an one front process (OFP). 
TFP is possible at some conditions (for example, at T rn < T m for a pure material) and do 
exist as a steady-state process in a definite region of control parameters. In this region there 
is one or a few OFP. For pure materials the TFP exists at — IA4I < A < A* (see Eqs. (19) 
and (23)), i.e. at T rn — (T m — T m ) < T < T m + (q — e)/c. When the initial temperature 
is T < T m < T m , the OFP of the transition SI — > S2 proceeds at the grain boundary 
51/52. When T m < Tq < T m , the OFP of the melting of 51 at the interface Sl/L is also 
possible. At To > T m the grain 52 can also be melted at the interface S2/L. The TFP 
can proceed faster than the corresponding OFP due to a heat transfer between melting and 
solidification fronts through a thin liquid layer. In order to initiate the TFP (especially at 
T m — (T m — T m ) < T < T m ) a liquid phase must be created inside the system. Then the 
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grain boundary S1/S2 splits into two interfaces, Sl/L and S2/L, and the TFP proceeds as 
a self- sustained process. 

It should be noted that the TFP can take place in a pure material in which there are 
several polymorphic modifications. Then a low temperature modification plays a role of 
the melted grain SI in Fig. 1 and its melting temperature T m is lower than the melting 
temperature T m for a high temperature modification, S2. In such a case the solutions given 
by Eqs. Q - are valid for A m = (T m — T m )c/q as a material parameter and (1 — aA m ) 
is replaced by another material parameter q/q, where q and q are the melting heats of the 
low and high temperature phases respectively. 

Generally, one can consider the TFP of the transition of a phase SI into another phase 
S2 through an intermediate phase L (which is not necessarily a liquid one) as an alternative 
to the OFP of the direct transition SI — > S2. 

The obtained Eqs. Q - together with Eqs. (J27|) - (J3U|) . which define relations 
between the Peclet numbers Pi, P2 and control parameters T , Co, give the continuous 
spectrum of solutions with the only free parameter (e.g., the velocity V). Therefore an 
additional equations, i.e. "selection" relationship, is needed in order to define the unique 
solution. This is similar to a well-known "selection" problem in dendritic growth [S|j. A 
search for the selection condition will be a subject of future investigation. The only point 
which might be stressed here is the following. 

The structure of the fronts is usually more complicated compared to pure parabolic ones 
due to possible cellular structures or due to finite size of the sample. In this case, two 
parabolas describe only a part of the moving fronts. The complete shape of the fronts must 
be subjected to additional boundary conditions. One can speculate on possible cellular 
structures with two moving interfaces which are shown in Fig. 3. These structures should 
appear due to the diffusional interaction between different cells and capillary effects play 
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also crucial role. The structures have central parabolic parts which are convex (Fig. 3(a,b)) 
or concave (Fig. 3(c)) and satisfy boundary conditions of zero heat- and mass- fluxes across 
boundaries of cells. The first structure (Fig. 3(a)) can be related to a melting processes and 
the other two ones, for example, to sintering of two solid grains 51 and 52 in liquid phase 
L. The structure in Fig. 3(b) may correspond to a sintering process with a supersaturated 
solid 51, while the structure in Fig. 3(c) with a concave central part may correspond to the 
case of undersaturated solid 51. 

In the process in the channel (or cell) the fronts velocity may be controlled either by the 
above mentioned cell boundary conditions or by the "selection" which takes place mostly in 
the parabolic region. Such a problem arises also in the selection of the growth velocity of 
the classical dendrite in the channel (the structure with one front) jH]. 

In this paper the solutions are obtained for two parabolic fronts in the two dimensions. 
Similar solutions can be easily found for two paraboloidal interfaces in the three dimensions. 

A few words about the solutions and (jHJ) for one convex and one concave parabolic 
interface might be outlined. The convex parabola and Eq. (P) are playing an important role 
in dendritic solidification (see, e.g., Ref. [S]) and must be important in "dendritic" melting. 
It can be supposed that concave parabola and Eq. (jHJ) are playing a role in such "doublon" 
structures with two convex parts and a concave part in between which are similar to the 
profile of 52/L-interface in Fig. 3(c). Whether or not those "doublon" structures exist in 
solidification or melting processes with one front? 

7 Conclusions 

A theoretical model for combined motion of two solid-liquid fronts with a liquid layer in 
between has been developed for pure materials and binary alloys. The temperature and 
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chemical compositions at the interfaces were supposed to be constants but different for 
fronts of melting and solidification. 

An exact solution which describes the steady-state motion of two co-focal parabolic in- 
terfaces has been found. It is shown that there exist two types of solutions with two convex 
and with two concave parabolas. The relations between the Peclet numbers and the control 
parameters of the process for both types of solutions are obtained. 

As usual in theories of the steady-state growth, the continuous spectrum of solutions 
with one free parameter exists. The unresolved problem of "selection" of unique solution is 
briefly discussed. 
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FIGURE CAPTIONS 

Figure 1: Two types of combined moving parabolic fronts: convex (left) and concave 
(right). 

Figure 2: The dependence of the reduced temperature A on the Peclet number P 1 for 
convex and concave fronts. The insert represents the region of small A. Parameters are: 
A m = 0.01 and a = 0.25. 

Figure 3: Three possible cellular structures of two combined moving interfaces with 
convex [(a) and (b)] and concave (c) central parabolic part. 51 and S2 are melting solid and 
growing solid, and L is the liquid phase. The direction of motion is defined by the arrow of 
V. Dashed regions in (a)-(c) show the central parabolic part of interfaces. 
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